Influence of boundary conditions on quantum escape 
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It has recently been established that quantum statistics can play a crucial role in quantum escape. 
Here we demonstrate that boundary conditions can be equally important - moreover, in certain cases, 
may lead to a complete suppression of the escape. Our results are exact and hold for arbitrarily 
many particles. 



I. INTRODUCTION AND SETUP 

The question of how particles escape from a partially 
confining region has long been at the center of many 
experimental and theoretical studies, leading to a num- 
ber of profound discoveries in mathematical physics (for 
a recent review see Ref. [1]). Most of the theoretical 
progress has been made in the context of classical me- 
chanics, where different dynamical behaviors - regular, 
chaotic, or mixed - lead to different escape laws; typically 
quantified by the decay of the survival (i.e. non-escape) 
probability. On the quantum mechanical side, there has 
been a surge of renewed interest in the escape proper- 
ties of few-particle systems [2-7] , in connection with re- 
cent advances in experimental control and manipulation 
of a small number of ultra-cold atoms [8, 9]. In particu- 
lar, particle-particle interactions and quantum statistics 
have been shown to significantly influence the escape. 
Here, we extend the theory to account for the influence 
of boundary conditions. 

We consider a system of N non-interacting particles in 
one spatial dimension which are initially (t — 0) confined 
within the interval I = (0, 1) with a hard impenetrable 
wall placed at the origin. A schematic plot of the setup is 
shown in Fig.l for N = 5 particles represented as Gaus- 
sian wave packets. For positive times, t > 0, the particles 





FIG. 1: Schematic plot of the initial configuration for N = 5 
particles represented as Gaussian wave packets. 



are allowed to explore the positive real line. The quan- 
tum state of the system is described by the N particle 
wave function ^(xi, . . . , xjy, t), which satisfies the time- 
dependent Schrddinger equation, 
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stant and particle masses to unity, i.e., h = m\ = 
. ..mjv = 1. The initial wave function is such that 
*(xi,.. .,x N ,t = 0) = ]£xj £ I for all j e [1, AT]. 

Under time evolution, the wave function diffuses out- 
side I. One way to measure this is through the survival 
probability pW(t), i.e. the probability that all N quan- 
tum particles are in I at times t > 0. In terms of the 
wave function, the survival probability is defined as 
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Quantum statistics obeyed by the N particles - 
whether bosonic or fermionic - can have a dramatic ef- 
fect on the time-decay of the survival probability. More 
specifically, it was observed in Ref. [2] that the sur- 
vival probability P^ 2 \t) asymptotically decays to zero 
as ~ t~ e for bosons and ~ t~ 10 for fermions. This sur- 
prisingly simple, yet generic, result was obtained by an 
asymptotic expansion of the appropriate quantum prop- 
agator for large t. This observation was later general- 
ized in Ref. [4] for N particles, showing that P^ N ^(t) 
decays like ~ t~ 3N and ~ t~ N ( 2N + 1 ) f or bosons and 
fermions respectively. The faster fermionic decay was 
attributed to the effective hard-core interaction among 
fermions causing anti-bunching and hindering the recon- 
struction of the initial state. In both works [2, 4] Dirichlct 
boundary conditions were imposed at the hard wall, i.e. 
^(a;!,...,^,*)!^ = for all j e [1,N\. 

When a particle is constrained to move within a lim- 
ited volume, Dirichlet (2?) or Neumann (Af) boundary 
conditions (BCs) are usually employed in order to reflect 
local probability preservation at x = 0. It is well known 
however that the general solution to this requirement is 
described by a Robin (JZ) BC [10], defined as a weighted 
combination of the T> and N BCs 



- V )t>(x 1 ,...,x N ,t) -0, Vj, (3) 



Hamiltonian operator and we have set the Planck's con- 



d 

dxi 



and controlled by the parameter r\ g R, such that 77 = 
corresponds to the and r\ = 00 to the T> BCs. The pa- 
rameter r) has the physical interpretation of a phase shift 
of the wave function on reflection with the wall at the 
origin. Therefore, the following natural question arises: 
What is the influence of boundary conditions on the de- 
cay of the survival probability? More specifically, how 
does the survival probability p( w ) (i) depend on 77 ? This 
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we address in the current setting by deriving and ana- 
lyzing the exact quantum propagator. It is worth noting 
that 1Z BCs are commonly used in Sturm-Liouville type 
problems which appear in many contexts in engineering 
and applied mathematics including for example in acous- 
tics and in convection-diffusion processes. 

We motivate this investigation by considering the sim- 
plest case of a single particle, N = 1. The dynamical 
evolution of the wave function can be described by the 
quantum propagator K v (x, x', t) as 



$(x,t)= / K n {x,x',t)V(x',0)dx'. 



(4) 



The propagator [11] is a generalized function that satis- 
fies the time-dependent Schrodinger equation 
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the initial condition 



J K, 1 (x,x / ,t)=0, 



lim K n (x 1 x , t) = S(x — x ) 
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i] \ K v (x,x',t) 



= 0, 



(7) 



the K BC at x = 

d_ 
^dx 

and that vanishes as x — > oo at negative imaginary times, 
i.e., lim KJx,x',-iT) = for T > 0. 

X—i-OO 

The cases of t] = (JV BC) and -q = oo (V BC) are 
particularly simple. Using the method of images we can 
immediately write down the full single-particle propaga- 
tor as 

K (x,x',t) = Kf [ee (x,x r ,t) + K flee (x, -x',t), (8) 
Kooix^x'ft) = K hcc (x, x',t) - K iicc (x, -x',t) , (9) 



where 
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(10) 



is the free-particle propagator, i.e. the propagator de- 
scribing the motion on the real line without a reflecting 
wall at the origin. 

The asymptotic decay of P^ (t) is entirely determined 
by the long-time expansion of the quantum propagator. 
Thus, 



and 



K (x,x',t) 



+ 0(t~ 3 ' 2 ) 
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+ 0{t- 5 / 2 ), (12) 



entailing 



TV BC : P (1) - t~ l , 
£>BC : P (1) -r 3 . 



(13) 



So for the single-particle case, the V BC leads to faster 
escape. 

When addressing the case of two particles, N = 2, the 
issue of quantum statistics comes into play. For noninter- 
acting particles, there are two ways of taking into account 
bosonic (b) or fermionic (f) statistics: either (i) by ap- 
propriately symmetrizing the initial state of the system, 

¥ h ^(x u x 2 ,o) = 4= [ihWMx*) ± MX2)MX1)] , 
v2 

(14) 

with "+" corresponding to bosons and "— " to fermions, 
and propagating it with K v (xi, x' x ,t)K v {x2, x' 2 , t), or (ii) 
by introducing an effective, appropriately symmetrized 
two-particle propagator 

K^( Xl ,X2,x[,x' 2 ,t) 

= —7= K v (xi,x / 1 ,t)K v (x 2 ,x / 2 ,t) 



± K ri (x2,x[,t)K n (xi,x' 2 ,t) 



(15) 



and applying it upon the product state il>i( x i)' l l > 2(x2)- 
Indeed, the two-particle state at t > can be written as 

*W(x 1 ,X 2 ,t) 

K r) (xi,x' 1 ,t)K v (x2,x 2 , 0) dx\dx' 2 

i 2 



= / K^ h/i \x 1 ,x 2 ,x l 1 ,x 2 ,t)ip 1 (x' 1 )tp2{x 2 )dx' 1 dx 2 , 
J i 2 



(16) 



where the first equality corresponds to (i) and the second 
to (ii) . Being interested in the asymptotics of the survival 
probability (and inspired by Ref. [2]), we adopt approach 
(ii) , as it allows us to focus only on the kernel, as opposed 
to the full integral. In other words, the long-time decay 
of p( 2 \t) is determined by Kq^(xi,X2,Xi, x' 2l t) in the 
limit t — > oo. For rj = and 77 = 00 one obtains 
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(17) 



It is interesting to note that the asymptotic decay of 
the survival probability of two fermions with the TV BC 
is governed by the same exponent as that of two bosons 
with the T> BC. This simple observation clearly illustrates 
that BCs imposed at the perfectly reflecting wall at the 
origin are as important for the particle escape as quantum 
statistics. It also suggests that by tuning the value of 
r\ one may expect to probe the continuum of possible 
decay rates, ranging between those for the T> and TV BCs. 
However, as we shall show in the following sections, this 
is not the case, as the decay exponent turns out to be a 
discontinuous function of r\. 
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II. SINGLE-PARTICLE CASE 

Unlike in the special cases of rj = (TV BC) and 
77 — oo (T> BC), the quantum propagator for an arbi- 
trary r] can not be straightforwardly obtained using the 
method of images; alternative techniques have to be em- 
ployed. We begin our study of the single-particle survival 
probability by deriving an exact closed-form expression 
for K rj (x,x' ,t) : defined through Eqs. (5-7). 



A. Exact propagator 

One way to construct the single-particle propagator 
K„[x, x' , t) is by using a complete orthonormal set of 
eigenstates {4>k} of the Hamiltonian H^ l > , 



Id 2 
2dx 2 



<f> k (x) = E k 4> k (x) 



(18) 



satisfying the 1Z BC at the origin, — 77) 4>k{0) — 0. As 
it will become clear from the following discussion (and 
as originally noted in Rcf.[12]), there is an important 
spectral difference between the case of 77 > and that of 
rj < 0. We analyze the two cases in succession. 

In the case of 77 > 0, the energy spectrum of 
is continuous. Indeed, the corresponding eigensystem is 
given by 



<t>k{x) 



\ 



7T 1 



n2 



k 

sinffea;) H — cos(kx) 



(19) 



with Ek = k 2 /2 and k > 0. The Hamiltonian eigen- 
states for the Af and T> BCs are respectively obtained 

- cos(fca;) and lim <fik(x) = 



as the limits lim <fik(x) 
77— >o 



f sin(fcx). 



rj— loo 



It can be verified by direct integration that the set of 
the eigenstates {</>&} is orthonormal, 



4>k(x)(j)k> (x) dx = 5(k- k') , 



(20) 



and complete, 

Mx)Mx') dk = S(x - x 1 ) , ?7>0. (21) 



Consequently, the propagator, expressed in the basis 
formed by {4>k}i reads 



K v (x, x , t) = (x\(- 



<j) k {x)<t>k{x')e *dk, f]>0 



(22) 



In view of Eq. (19), the integral in Eq. (22) can be rewrit- 
ten as 

1 f + °° 1 - 2 

— / [cos(x + k) + cos(x^k)]e~2 lk f dk 



7] f + °°Ticos(x + k) - ksin(x + k) _i lk 2 t 



7T 

ll -12 



k 2 + rj 2 



e~2 t,c l dk 



with x± — x ± x'. Then, it is easy to see that the 
first term is nothing but the quantum propagator for 
the J\f BC, i.e., T\ = K (x, x', t). The correspond- 
ing evaluation of the integral in the second term yields 

I 2 = rjerfc (j^^j e nx + + ^ ltTl2 , where "erfc" denotes 
the complementary error function. Combining the two 
terms we obtain our first main result 



K^lx^x' , t) = Kq(x, x , t) 

x + x' + irjt 



•x] erfc 



'2it 



(23) 



In the case of 77 < 0, the set of eigenstates {4>k} no 
longer forms a complete basis [12]. Indeed, in addition 
to the unbound eigenstates, forming a continuous part of 
the energy spectrum, there exists an isolated bound state 



X{x) 



-\f)\x 



V<0; 



(24) 



i.e. 



corresponding to the negative energy — 7y 2 /2, 
^ 1%) = ~ \\x)- The state is normalized to unity, 
(x\x) = 1) an( l is orthogonal to the eigenstates of the 
continuous part of the energy spectrum, i.e., (x\4>k) = 
for any k. The wave function x( x ) i g localized in the 
vicinity of the hard wall, indicating that the BCs here 
correspond to an effective zero-range attractive force. 

The presence of the bound state modifies the complete- 
ness relation (cf. Eq. (21)) 



(j) k (x)(j>k(x / ) dk + x(x)x(x') = S(x - x 1 ) , 77 < , 

(25) 

so that the propagator expanded in terms of the complete 
basis ({<M,x) reads (cf. Eq. (22)) 



K v (x,x , ,t)= cj>k{x)Mx')e~- 2lk l dk 
Jo 

+ x(x)x(x')e^ 2t , 77 <0. (26) 

The evaluation of the last integral for 77 < proceeds 
in the same manner as in the case of 77 > 0, and the 
resulting expression for the full propagator coincides with 
that given by Eq. (23). This concludes our derivation of 
the single-particle propagator. 

At this point we remark that the problem of a single 
particle restricted to the positive real line by a hard wall, 
is equivalent to that of two hard-core particles on the full 
real line through a change of variables to center of mass 
and relative coordinates. Moreover, the BCs imposed in 
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FIG. 2: (Color online) Snapshots of the evolution of the density | "3/ (a;, £) | 2 of a Gaussian wave packet initially centered at 
q = 0.6 with width a = 0.1 at times t = t w , 2t w ,4t w , for the top three panels and at t = t a ,3t a , 6t a , for the bottom three. The 
values t w = qa/3 and t a = q/n are respectively the approximate times when the wave packets start interacting with the wall, 
and when the asymptotic regime sets in (see discussion in main text for details). The four curves in each panel correspond to 
the evolution under different propagators K v with 77 = 00, 0, 2, and —2 as indicated in the last panel. 



the former correspond to the quantum statistics obeyed 
by the latter allowing for the possibility of anyons [12]. 
The equivalence in the case of many particles is however 
more subtle with essential differences in higher dimen- 
sions and will not be discussed in the current paper. 



B. Short-time dynamics 

Knowledge of the exact propagator (23) allows us to 
study the full dynamical evolution of the single particle 
wave function ^(x,t). In particular, in order to under- 
stand the influence of BCs on the escape, it is instructive 
to apply K n to a spatially localized initial state given by 
a Gaussian wave packet 



*(x,0) 



(27) 



centered at q £ (3er, 1 — 3er) with spatial width a -C 1, 
such that *(x £1,0) 0. 

Fig. 2 shows snapshots at different times of the proba- 
bility density function \^f(x, t)\ 2 under time evolution by 
Krj for four different values of 77. We denote these four 
cases by: (i) r\ = 00, (ii) 77 = 0, (Hi) 77 = 2 and (iv) 
7/ = —2. For the first two simple cases (i) and (H) (corre- 
sponding to the blue and red curves in Fig. 2), the density 



can be written down explicitly as 

2 e 1+- 2 



O^fll y/l + T 2 



cosh =F cos t. 

l + T 2 V 1 + r 2 



(28) 



where we have used the scaled variables r = t/a 2 > 0, 
£ = x/cr, and 9 = q/cr, for the sake of brevity. The 
"— " and "+" signs in (28) correspond to cases (i) and 
(ii) respectively. Cases (Hi) and (iv) in the figure were 
calculated through numerical integration of Eq.(4). In- 
specting Eq.(28) we can identify the first exponential as 
the main time envelope and its denominator as being 
indicative of the diffusive process. The cosine term de- 
scribes the oscillations due to interaction with the wall 
and are clearly seen in Figure 2. Note that the amplitude 
of these oscillations is suppressed for small times. 

We now qualitatively describe the dynamics observed 
in Fig. 2 with the aid of Eq.(28) and identify three dis- 
tinctly different time scales. For times < t < t w , all 
four curves coincide as they diffuse freely until the wave 
packets start interacting with the wall. This occurs when 



3<tVT 



q, and so t w w gcr/3 for a <C 1. For times 



t > t w , oscillations appear due to reflection off the wall 
and interference occurs. The wall interaction however is 
different in each case due to the different BCs prescribed 
by 77. It is clear from the figure that (i) and (ii) are com- 
pletely out of phase, while (Hi) and (iv) closely follow 
curve (ii) and are only slightly out of phase. The oscil- 
lations spread until t ss t a , when there is only a single 
maximum of |\l/(a;,t)| 2 left in / for case (i) and a single 
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minimum for case (it). The time t a can be extracted from 
the cosine term in Eq.(28) by requiring its argument to be 
2n at x = 1. We thus have that t a « q/n for a -C 1. In- 
deed, at times greater than t„, the asymptotic regime sets 
in and there are no more oscillations. Cases (i) to (Hi) be- 
have essentially the same, as they diffuse at different rates 
- while case (iv) seems to "stick" to the wall rather than 
diffuse away. The difference between (in) and (iv) can be 
more clearly seen in the two GIF animation files included 
as supplementary material. We postpone further discus- 
sion on the long time behavior of the four cases for the 
next subsection and turn to the survival probability func- 
tion. Fig. 3 shows the survival probability (2) for the four 




FIG. 3: (Color online) The single particle survival proba- 
bility in I for the four cases described in Figure 2 using the 
same values for q = 0.6 and a — 0.1. The asymptotic regime 
sets in around t ~ t a — q/n w 0.191. 



cases considered above for times up to 6t n , with the same 
parameters used in Fig. 2. Cases (?) and (ii) were cal- 
culated analytically, while (Hi) and (iv) by numerically 
integrating Eqs.(4) and (2). An additional time scale to 
the three described above now becomes relevant. For 
times t < td the Gaussian wave packets have not spread 
outside I and so there is not decay and PW(t < t^) = 1. 
Hence we have that = (1 — q)/3<r. Note that tg can be 
smaller, larger or equal to t w , depending on the value of 
q. From Fig. 3 we notice that pW(t) is decreasing non- 
monotonically and shows strong oscillations. These are 
due to the spreading of the interference oscillations ob- 
served in Fig. 2. Moreover, there are a finite number of 
oscillations as for i > 1, the frequency in the cosine term 
of Eq.(28) goes to zero. We can estimate the number of 
oscillations by f, ta -4wdt — — — 1 for a <C 1. 

For times t ^> a 2 , the two exponentials and the hy- 
perbolic cosine of Eq.(28) are approximately equal to 1. 
Using this approximation and integrating over x G I we 
obtain an approximation for P^ (t) for cases (i) and (ii) 



P«(t) 



2cr 



a 

1 T 7T sm 
2q 



(29) 



Expanding for large times produces the expected asymp- 
totics described by Eq.(13). Note that for case (i) cor- 
responding to T> BC, the order i _1 term cancels exactly, 
so P( x )(t)~t- 3 . 



C. Asymptotic analysis 

In the limits 77 — ^ and 77 — > 00 we obtain the propa- 
gators with TV and V BCs respectively, as expected, with 
next-to leading order corrections: 



K v (x, x , t) — Kq(x, x , t) 



erfc 



-i(x + x 1 ) 



2t 



v + o( v 2 ), 



(30) 



and 

K n (x,x',t) = K 00 (x,x',t) 
2i(x + x') 



Tft 



K fiee (x,-x',t) + 0( V - 2 ). 



(31) 



In the limit of rj — > —00 however we obtain 

K v (x,x',t) = -2|r;|e- | " l(a;+x ' )+ ^ t '' 2 
+ K 00 (x,x',t) + 0(r 1 ~ 1 ), 



(32) 



which implies that for negative (but finite) values of 77, 
the single particle survival probability saturates at some 
constant, rather than decaying to zero, as observed in 
Fig. 3. In the limit this constant vanishes and = 

To understand this better we expand Eq.(23) for large 
times t to get 



K n (x,x',t) = (M - 71 )e^ x+x ')+^ 2 
_ (1 + i)(l + xrj)(\ + x'ri) 

y^ V 2 t 3 / 2 



+ o(r 5 / 2 ). 



(33) 



Note that the first term vanishes when 77 > 0. Eq.(33) 
confirms that for 77 > the asymptotic single particle 
survival probability decays like ~ t~ 3 (as for T> BC) and 
for j] < the decay saturates at a constant C. We con- 
centrate on the first term in Eq.(33) and apply it to a 
Gaussian initial state as in (27) to get that 



C = A\r]\ae~ 2q ^ ( 1 



-2\V\ 



O(o 2 



(34) 



For 7/ = —2 and a — 0.1 we have that C w 0.126 which 
is in very good agreement with the asymptotic decay of 
case (iv) in Figure 3. 

We remark that Eq.(33) can be seemingly misleading, 
since for 7; = (corresponding to TV BC) the second term 
diverges like ~ —^/2i/(y/ : Krj 2 t i / 2 ). Indeed this implies 
that the limit 77 — > does not commute with that of t — > 
00. In fact we notice that Eq.(33) reproduces the correct 
asymptotics prescribed by Eq.(ll), only if 77 scales like 
\f\Ji. This can be understood in the following way: for 
7] < 0, the bound state x becomes completely delocalized 
in the limit 77 — > 0~. One can then imagine that \ acts 
like an "?7-thin carpet" hindering the particle diffusion, 
thus giving a slower decay of P*- 1 ' ~ 

In brief, the exponent a of the asymptotic power law 
decay PW - /T Q experiences discontinuous jumps from 
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for 77 < 0, to 1 for r\ = 0, and 3 for r\ > 0. This unex- 
pected discontinuity has interesting repercussions when 
considering the escape of N > 2 particles from the inter- 
val / = (0, 1). These are discussed in the next section. 



III. MANY-PARTICLE CASE 

In the case of two particles, the initial state is given by 
Eq.(14) with 



1 



(35) 



and the wave function vp( b / f ) (xi,X2,t) at a later time 
t > is determined by Eqs.(15) and (16). We have pre- 
viously discussed the asymptotic behavior of P^ 2 ' for Af 
and V BCs (see Eq.(17)) and now turn to the general 
case of arbitrary 77. This is achieved by expanding the 
propagator (15) in the limit of large t: 



Bosons, 77 < 
Fermions, 7/ < 

Bosons, 77 > 
Fermions, 77 > 



K ( f ) ■ 

K (b) 



' const 

t -3/2 



-3 



> PW 
p( 2 ) ^ 
p( 2 ) ~ 



/.-6 



const . 

3 



(36) 



It is interesting to note that the anti-symmetry of the 
fermionic wave function does not allow for a bound state 
and thus P^ decays to zero. This completes the picture 
for the two particle case. 



TABLE I: The asymptotic survival probability for one, two, 
and N particles, as well as different BCs. 



K BC 


1 


2 bos 


2 fer 


N bos 


N fer 


77 < 


const 


const 


t~ 3 


const 


t (l-N)(2N~l) 


T] = 


t- 1 


t~ 2 


t~ 6 


t -N 


.-JV(2JV-1) 


77 > 


t- 3 


t- 6 


t- w 


t -3N 


^-JV(2JV+1) 


Tj = ±00 


t~ 3 


t~ 6 


t -io 


t -3N 





The propagator (15) can be generalized for TV > 2 
particles in a standard way [13]: 

K^ t \x l ,...,x Nl x' 11 ...,x l N ,t) = 
1 N 

~7W] ^ t e ii,..-,iN-^ri{ x l,X il ^t)...K ri (xN,X iNl t), 
il,...,ijV=l 

(37) 

where £i 1) ...,i N is the totally antisymmetric (permutation) 



symbol in the case of fermions and is identically 1 for 
bosons. 

As done before, the asymptotic decay is obtained by 
expanding Eq.(37) for large times t. We present the re- 
sults for N = 1,2, and for the general case in Table I al- 
lowing for a direct comparison. In particular, we observe 
that the case of 77 > coincides with that of 77 = ±00 
corresponding to T> BC. For the case of 77 < 0, the proba- 
bility decays to a constant for bosons due to the presence 
of a bound state. Significantly, we observe that the case 
of 77 = 0, corresponding to JV BC, is special as it sepa- 
rates the other cases in a discontinuous fashion for both 
bosons and fermions. We also observe that fermions al- 
ways escape faster than bosons. As we have previously 
pointed out, the asymptotic decay for two fermions with 
Af BC is the same as for two bosons with T> BC. However, 
this correspondence no longer holds for N > 2. Finally, 
we note that the survival probability of N fermions with 
77 > (or 77 = ±00) has the same decay exponent as that 
of N — 1 fermions with 77 < 0. 



IV. CONCLUSIONS 

We have addressed the influence of boundary condi- 
tions on the escape of N indistinguishable particles in a 
one-dimensional setting, and have shown that they are 
equally important as the prescribed quantum statistics. 
To this end, we have derived an exact closed form expres- 
sion for the single-particle propagator on the positive real 
line in the presence of Robin boundary conditions with a 
single control parameter 77. This expression generalized 
existing results restricted to Dirichlct boundary condition 
(7/ = 00) and unveiled new non-trivial scenarios where es- 
cape may be completely suppressed. Moreover, we have 
found that the exponent of the asymptotic power law de- 
cay of the survival probability is a discontinuous function 
of 77. Our results hold for an arbitrary number of particles 
and are summarized in Table I. 

In the light of recent atom-optics experiments ad- 
dressing the dynamics of a small number of bosons and 
fermions in one dimension [9], our findings may lead to 
new applications in the area of quantum control and state 
manipulation. On the theoretical side, it is also interest- 
ing to explore how higher dimensional waveguide- like [14] 
or network geometries [15], external potentials, particle- 
particle interactions may affect the dynamics and in par- 
ticular quantum escape. 
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